Preprint typeset in JHEP style - PAPER VERSION 



DPNU-04-24 
IU-MSTP/64 
hep-lat/0412042 



A local formulation of lattice Wess-Zumino model 
with exact U{1)r symmetry 



Yoshio KIkukawa 

Department of Physics, Nagoya University, Nagoya ^6^-8602, Japan 
E-mail: kikukawaOeken . phys . nagoya-u .ac.jp 

Hiroshi Suzuki 

Institute of Applied Beam Science, Ibaraki University, Mito 310-8512, Japan 
E-mail: hsuzukiOmx . ibaraki .ac.jp 

Abstract: A lattice Wess-Zumino model is formulated on the basis of Ginsparg- Wilson 
fermions. In perturbation theory, our formulation is equivalent to the formulation by 
Fujikawa and Ishibashi and by Fujikawa. Our formulation is, however, free from a singular 
nature of the latter formulation due to an additional auxiliary chiral supermultiplet on a 
lattice. The model posssesses an exact U(l)/j symmetry as a super symmetric counterpart 
of the Liischer lattice chiral U(l) symmetry. A restration of the supersymmetric Ward- 
Takahashi identity in the continuum limit is analyzed in renormalized perturbation theory. 
In the one-loop level, a supersymmetric continuum limit is ensured by suitably adjusting 
a coefficient of a single local term F*F. The non-renormalization theorem holds to this 
order of perturbation theory. In higher orders, on the other hand, coefficents of local terms 
with dimension < 4 that are consistent with the U(l)ij symmetry have to be adjusted for 
a supersymmetric continuum limit. The origin of this complexicity in higher-order loops 
is clarified on the basis of the Reisz power counting theorem. Therefore, from a view point 
of super symmetry, the present formulation is not quite better than a lattice Wess-Zumino 
model formulated by using Wilson fermions, although a number of coefficients which require 
adjustment is much less due to the exact U(l)ij symmetry. We also comment on an exact 
non-linear fermionic symmetry which corresponds to the one studied by Bonini and Feo; 
an existence of this exact symmetry itself does not imply a restoration of supersymmetry 
in the continuum limit without any adjustment of parameters. 
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1. Introduction 



There has been a renewed interest on non-perturbative formulation of supersymmetric 
theories via a spacctime lattice [l]-[9] in these several years [10]-[21] (for a recent review 
with a complete list of references, see rcf. [22]). One major idea in these recent studies is 
to keep a part of the supersymmetry algebra manifest and infer that this exact symmetry 
is strong enough to ensure a fully supersymmetric continuum limit without any (or with 
a small number of) adjustment of parameters [10, 15, 16, 19]. This general strategy, 
which is also common to some of past attempts [2], has achieved fair success, typically 
for lower dimensional supersymmetric theories with an extended supersymmetry (besides 
a potential problem of positivity of the measure). An extended supersymmetry allows a 
sub-algebra which is consistent with a lattice construction and, due to the lowcr-ness of 
dimensionality, the number of relevant operators, which potentially violate supersymmetry 
in the continuum limit, is small. So often a supersymmetric continuum limit is achieved 
without any adjustment of parameters. 

Another kind of approaches is to abandon a manifest supersymmetry of a lattice model 
from the onset and achieve a supersymmetric continuum limit by adjusting parameters in 
the model. This is an approach advocated in refs. [3, 4, 5, 7] and this has been, in our 
opinion, only a realistic approach to date for A?^ = 1 supersymmetric theories in 4 dimen- 
sions. Here, again, some exact global symmetries on a lattice can be useful [8, 9] to reduce 
the number of parameters which require adjustment. For numerical simulations along this 
kind of approaches, sec ref. [23].^ 

In this paper, we adopt the latter attitude and study a lattice formulation of the 
4 dimensional A^ = 1 supersymmetric Wess-Zumino model [24]. The Wess-Zumino model 
is asymptotic non-free and thus the continuum limit, as a fundamental theory, is expected 
to be trivial. Nevertheless, it is meaningful to consider the model as an effective theory 
with an ultraviolet cutoff. In a sense, this model is more difficult to formulate on a lattice 
than supersymmetric Yang-Mills theories because a quadratic divergence in mass terms of 
scalar fields is expected to be prohibited only with presence of an exact supersymmetry. It 
is thus deserve to study in its own right. We formulate the model on the basis of Ginsparg- 
Wilson fermions [25]-[29]. This kind of formulation has been pursued by Fujikawa and 
Ishibashi and by Fujikawa [11, 12, 14]. In our notation, their formulation is expressed as 

^ = - laD2r'Dix - l(l^*D2cl) + F*{l - ^aD2)-'F 

X 

+ -mx^CP+x + -m*x^CP-X + 'mFcj) + m*F*(f)* 

+gX^Ccj>P+x + gFcP^ + 9*x^C<j>*P^x + 9*F*<f>*^], (1.1) 

where (%, ^, F) denote the chiral multiplet of the Wess-Zumino model and Di and D2 are 
lattice difference operators which will be defined below. In ref. [11], explicit perturbative 

^ There exist more ambitious approaches which aim an exact full supersymmetry on a lattice [6, 13, 20]. 

Another interesting observation is that a supersymmetric continuum theory is automatically restored if a 
convergence behavior of Feynman integrals in a lattice model is moderate enough [14, 17]. 
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calculations in one-loop order were carried out and it was found that, in the one-loop level, 
effects of supersymmetry breaking in the model appear only in wave function renormaliza- 
tion factors of the chiral multiplet, thus the violation of supersymmetry is rather moderate 
in the one-loop level. 

One can carry out perturbative calculations on the basis of the action (1.1) without 
any problem. However, as pointed out in ref. [12] (sec also ref. [30]), the action itself 
is singular because the operator 1 — {\/2)aD2 always has zero modes. This also implies 
that the kinetic operator in eq. (1.1) is non-local. Thus the meaning of the model in a 
non-perturbative level is not clear. 

In this paper, we first formulate a non-singular local lattice action for the Wess- 
Zumino model which is, in perturbation theory, equivalent to the formulation based on the 
action (1.1). This is achieved by introducing a non-dynamical auxiliary chiral multiplet 
on a lattice which decouples in the continuum limit. Due to the Ginsparg- Wilson relation, 
when m = 0, our model possesses a lattice analogue of the U(1)_r symmetry which is 
supersymmetric counterpart of the Liischer lattice chiral U(l) symmetry [29].-^ This is the 
symmetry that was pointed out for a free theory in ref. [9]. These are contents of section 2. 

Next, wc study a restoration of supersymmetry in the continuum limit by using a 
lattice version of the Ward-Takahashi identity. We carry out an explicit one-loop eval- 
uation of the supersymmetry breaking term in the Ward-Takahashi identity and observe 
that effects of supersymmetry breaking in the present model appear only in the wave 
function renormalization of the auxiliary field F in the continuum limit (in the one-loop 
level). We then present a general argument for higher loop contributions of supersymmetry 
breaking in renormalized perturbation theory. Unfortunately, the general power counting 
argument which is based on the Reisz theorem [31, 32] indicates that all supersymmetric 
non-invariant local terms with the mass dimension < 4 are radiatively induced, unless a 
term is forbidden by the U(1)r global symmetry that is manifest in our formulation. We 
clarify the reason why the one-loop result is so simple and higher loop corrections are 
expected to be destructive. In terms of the power counting theorem, the supersymmetry 
breaking term, which is a consequence of a violation of the Leibniz rule on the lattice, 
behaves as a non-derivative coupling in one-loop diagrams while it behaves as a derivative 
coupling in higher loop diagrams. This peculiar behavior of the supersymmetry breaking 
term makes the situation in higher loop diagrams involved. As a conclusion, from a view of 
supersymmetry restoration, our formulation is not quite better than the formulation based 
on the Wilson fermion [3] , although some of super non- invariant local terms are prohibited 
by the exact U(1)r symmetry. (Sec. 3) 

In the final part of this paper, we will comment on an exact non-linear fermionic sym- 
metry in our formulation which corresponds to the symmetry recently studied in ref. [21] in 
the context of the Fujikawa-Ishibashi formulation. This symmetry is nothing but the "lat- 
tice supersymmetry" utilized in ref. [7] for 2 dimensional Wess-Zumino model. As noted in 
ref. [7] and as indicated from the results of ref. [11] and of ours, a presence of this symmetry 
itself does not imply an automatic restoration of supersymmetry in the continuum limit 

^In the action (1.1), this U(l)ij symmetry is trivially realized as SaX = +*ck75X) = — 2ia0 and SaF = 
+4AaF. See section 2. 
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without any adjustment of parameters (although the non-hnear symmetry reduces to the 
standard supersymmetry in the classical continuum limit). We clarify this point. 
Throughout this paper, the lattice spacing will be denoted by a. 

2. The model 

2.1 Action and its symmetries 

Our staring point is the chiral invariant lattice Yukawa model of rcf. [29]: 

+2g* + *) (<^* + ^) P_ (V' + ^-j } , (2.1) 

where the field is a non-dynamical auxiliary fermionic field and P± = (1 it 75)/2. In this 
expression, we have shifted the scalar field as, cf) ^ (j) + m,/{2g) to generate mass terms for 
fermions. As the lattice Dirac operator D, we adopt the overlap-Dirac operator [26, 27] 
defined by^ 

D = i{l-A(At^)-V2}, A = l-aD^, D^ = \{^^{dl + d^)-adld^}, (2.2) 

which obeys the Ginsparg- Wilson relation 75D + D75 = aD^^D [25]. Thanks to this 
relation, the action with m = is invariant under a lattice chiral transformation of the 
following form [29] 

5atp = iaipil - ^aD)-f5 + ia^js, Jq* = ia'ijj^aDj5, 

^a4^ — — 2ia(/>, 5a,(t>* = 2ia(f)*, (2-3) 

where a is an infinitesimal real parameter. This transformation is designed so that a sum of 
fields, say + transforms in a standard way, 6aitp + ^) = iaj^iip + ^)- Thus a breaking 
of this chiral symmetry due to the presence of mass terms has a simple structure as in the 
continuum theory. The auxiliary fermion ^ is introduced to make a chiral transformation 
of this standard form and the Ginsparg- Wilson relation (which implies a non-standard 
chiral property of the lattice Dirac operator D) compatible. 

To define the Wess-Zumino model, we need to reduce degrees of freedom of the Dirac 
field ip to the Majorana one. Since the chiral projectors P± = (1 it 75)/2 in the Yukawa 
interaction term are ordinary ones, the Majorana reduction (see ref. [11]) can be applied 
straightforwardly. Namely, we make substitutions^ 

= (x + iv)/^, ? = ix^C - iv^C)/V2, 
^ = {X + iY)/ V2, = {X^C -iY^C)/V2, (2.4) 

^dfj,f{x) = {f{x + aft) — f{x)}/a and d*f{x) = {f{x) — f{x — ajl)}/a are the forward and backward 
difference operators, respectively. 

'^C is the charge conjugation matrix which satisfies C7jaC~^ = — 7J, C75C~^ = 7J, C'^C = 1 and = 
-C. 
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in the action. Noting relations 

{CDf = -CD, {CP±) = -CP±, (2.5) 

we find that the action decomposes into two independent systems. By taking only terms 
including x and X, we have^ 

+g*[x'' + X^)c(<t>* + ^)P-{X + X)]. (2.6) 

When TO = 0, the action is still invariant under the chiral transformation 

5aX = ^"75(1 - ^a-D)x + ia-y^X, 5aX = ia-f5^aDx, 

Sacp = -2ia(l), Saf = 2ia(j)*. (2.7) 

Eq. (2.6) provides a part of our lattice Wess-Zumino model, kinetic terms of fermions and 
Yukawa interaction terms. 

We next introduce bosonic superpartners of fermion fields, {(p, F) and (<I>, JT), and seek 
an appropriate free action which is invariant under a certain "lattice supersymmetry" . As 
the form of this fermionic transformation, we postulate 

S,x = -V2P+{D^cl) + F)e - V2P_{Dicl>* + F*)e, 

5, ct> = V2e^CP+x, = V2e^CP.x, 

6eF = V2€^CDiP+x, SeF* = V2e^CDiP-x, (2-8) 

and 

6, X = -V2P+{Di^ + J^)e - V2P_{Di^* + J^*)e, 
(5,$ = V2e^CP+X, 6e^* = V2e^CP-X, 

6eJ^ = V2e^CDiP+X, d,J^* = V2e^CDiP_X. (2.9) 

In this expression, e is a 4 component Grassmann parameter and we have used a decom- 
position of the Dirac operator, D = Di + D2, where 

Di = ll.id; + d,){A^A)-y^ = i{l - (1 + \aXd,){A^A)-'/'}. (2.10) 

Note that with respect to spinor space, Di and D2 have different structures. In particular, 
we have {75,1)1} = and [75,1)2] = 0. In terms of this decomposition, the Ginsparg- 
Wilson relation is expressed as 

2D2 = ai-Dl + Dl), (2.11) 



^This Majorana reduction, in a level of the functional integral, corresponds to the prescription of ref. [33]. 
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and as a consequence, we have relations 

75(1 - ^aD)js{l - ^aD) = 1 - ^aL>2, 75(1 - ^aD)j5D = Di, (2.12) 

which win frequently be used below. It is also understood that the 4x4 identity matrix 
in operators Df and D2 is omitted when these operators are acting on bosonic fields. It 
is then straightforward to see that the following free action is invariant under eqs. (2.8) 
and (2.9): 

X 

--X^CX - -(J^$ + J^*<^*) 
a a 

+^mx^CP+x + ^m*x^CP_x + mF^ + m*F*i>*], (2.13) 

where we have introduced abbreviations 

X = X + X, 4> = ^ + ^, F = F + r. (2.14) 

The combinations {x,(f),F) and {X,^,J^) are regarded as chiral multiplet in the lattice 
model. In particular, we refer (X, ^,J^) to as the auxiliary chiral multiplet which is char- 
acteristic in the present lattice formulation. 

We note that the free action with m = possesses three types of U(l) symme- 
try [9]. The first is a rather trivial one acting only on bosonic fields and is defined by the 
transformation: 

SaX = 0, 5aX = 0, 

SaF = —iaF, botT = —iaJ^, (2-15) 

where a is an infinitesimal real parameter. This remains the symmetry of So even for m 7^ 0. 
The second one is nothing but the Liischer chiral symmetry, (2.7) with 5a4> = 0, = 0, 
5aF = and = 0. Thirdly, somewhat surprisingly, the bosonic sector of Sq with m = 
possesses an analogous U(l) symmetry to eq. (2.7): 

SaX = 0, 5aX = 0, 

Sa<P = +ia{{l - ^aD2)(l) - 2"-^*} + ^<^*' = +ia{-aD2(t) + 2"-^*}' 

6aF = +ia{{l - ^aD2)F - ^aDfcj)*} + iaJ^, b^T = +ia{]^aD2F + ^aD?0*}, 

(2.16) 

due to the Ginsparg- Wilson relation. The lattice action Sq is not invariant under a uniform 
rotation of the complex phase of bosonic fields, (/>, F, <I> and JT, due to the presence of 
terms FD2(t) and F*D2(t)* ■ The above provides a lattice counterpart of this uniform phase 
rotation of bosonic fields under which the free action Sq with m = is invariant. Using a 
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linear combination of the above three U(l) symmetries, it is possible to define the U(l)ij 
symmetry [9] in the interacting system, as we will see below. It is worthwhile to note that 
a sum of transformations (2.15) and (2.16) takes the following simple form when acting on 
tilded variables: 

SaX = +^a75X) Sa4> = +ia4), 5aF = +iaF. (2-17) 
Next we postulate a form of the interaction term as 

Sint. = a''J2{sx''C4>P+X + gF^^ + g*X^C4>*P-X + 9*F*r^], (2.18) 

where we have defined interaction terms by taking tilded variables (2.14) as unit, because in 
this way we can relate our formulation to the Fujikawa-Ishibashi formulation. This way of 
construction of interaction terms is also a natural generalization of the Yukawa interaction 
in eq. (2.1). We then find that the full action S = Sq + S'int. is "almost" invariant under 
the fermionic transformations (2.8) and (2.9). In fact, after some algebra using the Fierz 
identity,^ we obtain 

S,S = -a''Y,x'^CV2[gP+{2^D^^-D^^^)e + g*P^{2^*D^^* -Di^*^)e] 

X 

= -a^Y^X^CALe. (2.19) 

X 

We emphasize that this breaking could vanish if the Leibniz rule was valid for the lattice 
difference operator D^. In summary, the lattice action for the Wess-Zumino model 

X 

a a 
+^mx^CP+x + ^'rn*x^CP-x + 'mF4> + m*F*4>* 

+gx'^C^P+X + gF4>'' + g*x^C^*P-X + g*F*4,*^Y (2.20) 

is invariant under the lattice super transformation (2.8) and (2.9) up to the breaking 
term (2.19). 

For the action (2.20), we can define two types of exact global "symmetries". The first 
is eq. (2.15) which yields on tilded variables 

5^X = 0, Sa4> = +ia4>, SaF = -iaP. (2.21) 

This is not a symmetry when g ^ 0, but may be regarded as a "symmetry" if we simulta- 
neously rotate the coupling constant according to 

dag = -lag. (2.22) 

®We use the identity (x^C'x)x = -(x^C'75X)75X- 
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Another is a lattice counterpart of the U(l)/j symmetry which is given by a Unear 
combination of the above three U(l) transformations: 

SaX = +^075(1 - ]j^a.D)x + ia-y^X, S^X = +ia'y5^aDx, 
5a4> = — 3ia^ + ia{(l — -aD2)4> — + ia$, 

Sa^ = -3ia^ + ia{^aD2^ + 

5aF = +2,iaF + ia{{l - ^aD2)F - ^aDfcp*} + iaF, 

5aJ^ = +SiaJ^ + ia{^aD2F + ^aDl4>*}. (2.23) 

On tilded variables, this U(l)ij transformation takes a simple form 

SaX = +io^l5X, Sa4> = -2ia4>, 6aF = +AiaF. (2.24) 

The action S with m = is invariant under this transformation and this may also be 
regarded as a "symmetry" even for m 7^ if we transform the mass parameter according 
to 

5am = —2iam. (2.25) 

Eqs. (2.23), (2.24) and (2.25) correspond to the U(l)ij transformation of the continuum 
Wess-Zumino model. In fact, if we define the U(l)ij transformation of the parameter of 
the super transformation as 

5ae = — 3ia75e, (2.26) 

then from eqs. (2.8) and (2.9) it can be confirmed that [(5e,(5Q] = holds on all field 
variables. The above two "symmetries" play an important role when we consider a structure 
of radiative corrections in the present model. 

In our non-singular local action (2.20) with interactions, the chiral U(l)ij symmetry is 
realized as an exact symmetry. The no-go theorem of ref . [30] on a chiral invariant Yukawa 
interaction of the Major ana fcrmion is evaded in our formulation due to an introduction of 
the auxiliary field (s). We further clarify this point in the next subsection. 

2.2 Perturbative equivalence to the Fujikawa-Ishibashi formulation 

In perturbation theory, the above system S is completely equivalent to a lattice Wess- 
Zumino model formulated in refs. [11, 12, 14], i.e., eq. (1.1). A formal way to see this 
equivalence is to rewrite the action S in favor of tilded variables and of the auxiliary 
multiplet {X,^,T): 

^ = - X^)CD{x -X) + (0* - ^*)Dl{4> - ^) 

X 

+(^* _ jr*)(F - JT) + (F - J^)D2{^ - $) + {F* - J^)D2i4>* - $*) 

--X^CX --(J^^ + J^*^*) 
a a 
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1 I _ _ 

+-mx^CP+x + ^'rn*x^CP-x + + m*F*4)* 

+gx^C4>P+x + 9F4>^ + g*x^C4>*P-x + g*F*^*^y (2.27) 

If we perform integrations over the auxiliary chiral multiplet, X, $ and we have the 
effective action for tilded variables: 

X 

+\mx-CP,x + \m*x'CP.x + mF~^ + 

+gx^C~4>P+x + 9F4>^ + g*x^C4>*P-x + g*F*4>*^y (2.28) 

This is, if we identify tilded variables as basic field variables, nothing but the action (1.1). 
Associated to the integration over the auxiliary chiral multiplet, we have 

Pi{c{D - -)} det-'\-{D2 - -)|, (2.29) 

where the first factor comes from an integration over the 4 component fermionic spinor X 
and the second comes from an integration over complex bosonic scalars $ and the 
operator in the latter factor therefore does not contain 4x4 identity matrix in spinor 
space. In a formal sense, these two factors are cancelled out, because the relation 

75(L>--)75(^--) = --(^2--) (2.30) 
a a a a 

holds (recall eq. (2.12)) and thus we have 

det^/D--) =det4( — (r»2--)), (2.31) 
I a) la a ) 

by noting the fact that the right hand side of eq. (2.30) contains the 4x4 identity matrix. 
Therefore, we see 

Pf{c7(D - ^)} = detV2 (^^gti/2|^ ~ -} = det^^^ C det\^^{D2 - ^)}. (2.32) 

This cancels the contribution from the bosonic fields. The system S is thus equivalent to S 
after integrating out the auxiliary chiral multiplet, X, $ and J^. 

This argument for the equivalence between 5* and S, however, is valid in a formal 
sense, because the operators D ~ ^ and thus D2 — | which appear in various places in 
the above expressions always have zero modes when the lattice volume is infinite [12, 30]. 
The kinetic operators of X and of have zero eigenmodes when tilded variables are 

kept fixed. Thus the integration over the former gives zero and the latter gives infinity. On 
the other hand, kinetic operators in the effective action S contain the factor (1 — ^aD2)'~^ 
which is a singular operator. 

To see what is really happening here, it is instructive to consider the case of m = 5 = 0. 
In this case, integrations over Grassmann variables yield 

JJdx(x)dX(x)e-"'^-{5X^CDx-ix^cx} ^ pfic-i^iPfj-^c}, (2.33) 



/: 
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which is not singular in any sense. On the other hand, if we perform the integration over X 
first while keeping x = X + we have instead 

pdc{D--)} /'j]dx(x)e-'^'S^{5X^f^a-|«^2)-ii?ix}^ (2.34) 

which is of a structure of x oo, although nothing is wrong with the whole integral (2.33). 
The above argument simply shows that we are observing a non-singular object in an unnec- 
essarily singular way. A similar argument is applied to the full action (2.20). Integrations 
over field variables do not produce any singularities. Only if we observe the integrations 
in a wrong way, seemingly singular natures as in eqs. (2.28) and (2.29) emerge. The full 
action S (2.20) and the action S (2.28) are thus inequivalent by a singular quantity.'' 

Nevertheless, we can infer that our formulation based on 5 and that based on S are 
equivalent in perturbation theory. The point is that the free propagators among tilded 
variables 



-Di + (1 - ^a£)2)(m*P+ + mP-) ^_4^ 
|L>2 + (1 - ^aD2)m*m " 



-D2 + (1 — ^aD2)rn*m 

{^{x)F{y)) = {F{x)~4>{y)) = 



(1 - ^aD2)m* _4 
ii^2 + (l-iaD2)— " 



(0*(x)F*(y)) = (F*(x)0*(y)) = '"^'^"^^^.^^ 



ID2 + (1 - ^aD2)m''m 
[1 - \aD2)lD2 
ID2 + (1 - \aD2)m*m^ 



{F{x)F*iy)) = {F*{x)F{y)) = ' ,\ a-%,y, (2.35) 



which are directly obtained by using S, are identical to those obtained by using S formally, 
i.e., by neglecting a singular nature of kinetic terms. Moreover, interaction vertices of S are 
identical to those of S (in fact we have constructed S so that this is the case). Therefore, 
perturbative calculations based on S and that based on S give rise to completely identical 
answers for correlation functions which have tilded variables in external lines. In a sense, 
our non-singular local lattice action 5 provides a natural justification for a prescription of 
refs. [11, 14] which utilizes the above form of propagators and interaction vertices of S. 
Of course, we think our formulation which includes the auxiliary chiral supermultiplet is 
superior at least formally, because it is manifestly free from singularities and it may have 
a meaning even as a non-perturbative formulation. 

3. Supersymmetric Ward-Takahashi identity and its breaking 
3.1 Derivation of a lattice Ward-Takahashi identity 

We consider a structure of radiative corrections with the present lattice formulation of 

the Wess-Zumino model. As noted in the preceding section, in perturbation theory, our 

'^In this way, the no-go theorem of ref. [30] for a chiral invariant Yukawa interaction of the Majorana 
fermion is evaded by introducing auxihary fields. 
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formulation is equivalent to the formulation of refs. [11, 12, 14]. One- loop radiative cor- 
rections in the latter formulation, in view of a realization of supersymmetry, had been 
extensively studied in ref. [11]. Here we study this issue in the continuum limit by using a 
Ward-Takahashi identity. 

For a systematic study, it is quite helpful to introduce the one-particle irreducible 
(IPI) effective action F. Following the standard procedure, we introduce external sources 
for each elementary fields 

^source = a'^J2^'^xX+J<t>(l^ + Jr4>* + JFF+JF*F* + JxX + J^^ + J^*^* + Jj.J^+Jj.*J^*}. 

X 

(3.1) 

We also introduce an external source K for a symmetry breaking of the action, SeS, and a 
source L for a symmetry variation of 5^3, 6^t5eS, and so on. Including these latter kind of 
external sources K, L. . . only, we define the total action 

Stot. = -S-a^^{K„(x^CAL)„ + L„,;35a(x^C'A% + M„;3,^5a(5;3(x^C'AL)^ + - • •}, (3.2) 

where 6a stands for the symmetry variation with the transformation spinor parameter is 
removed: 

Se = CaSa- (3.3) 

The generating functional W of connected Green's functions is then defined by the func- 
tional integral 

f Y[ dx{x)d(f>{x)d(f)*{x)dF{x)dF%x)dX{x)d^x)d<^*{x)dT{x)dr*{x) 

X 

'S'tot. 'S'source 

We then apply the Legendre transformation to W and change independent variables from 
external sources for elementary fields {J^, • • •) to the corresponding expectation values of 
elementary fields ((x)) " " ")■ what follows, we denote expectation values by their original 
name as {x) X <i^(l so on for notational simplicity. Wc do not apply the Legendre 
transformation with respect to the sources {K, L, M, . . .). In this way, we have the IPI 
effective action 

r = r[x, 4>, ^*,F, F*,X, J^, J^*;K, L, M, . . .], (3.5) 

which is a generating functional of IPI Green's functions which include additional vertices 
coming from the second term of eq. (3.2). 

Now, the action S is not invariant under the lattice super transformation (2.8) and (2.9), 
but it leaves the breaking (2.19). From this fact, we have the identity 



'^{J^SeX + J<l>5e<t> + J^*Se(t>* + JpSeF + Jp^SeF* 



X 

J2{KaeMx^CAL)a + L„,^e^5^,5«(x^CAL)^ + . . .}\ = o. (3.6) 
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This identity, in terms of the IPI effective action F, is expressed as 



- r^{V2P+{Di<t> + F)e + V2P_{Di(j)* + F*)e} 



dx 



X X 

- F-^{V2P+{Di<^ + J-)e + V2P_(Z)i$* + J^*)e} 

X 

X 

d 



a; 



a; 



that is a hnear equation of F. This is referred to as the lattice Ward-Takahashi identity. 
If the last two hnes were not present in this expression, the above equation simply states 
that the effective action is invariant under the lattice analogue of super transformation, 
eqs. (2.8) and (2.9). Thus, contributions of these lines, especially contributions from the 
term, 

that is the breaking of the supersymmetric Ward-Takahashi identity, will play a central 
role in our analysis below. Explicitly, this term is given by IPI diagrams with insertions of 
the operator*^ 

X 

= -a^ l]x^CV2{5P+(2^Di^-L>i^2)e + 5*P_ (2^*1^1^* -Di^*2^e}. (3.9) 

X 

The a — > limit of these IPI diagrams will be expressed by a local polynomial of field 
variables, because in this limit, the effect of our particular choice of regularization (the 
lattice artifact) should affect only local terms in the effective action F. Moreover, the op- 
erator (3.9) vanishes in the classical continuum limit (because the Leibniz rule holds in this 



^Generally, F contains IPI diagrams with multiple insertions of this operator. We will be interested 
in, however, terms of F that are linear in the external source K and consider IPI diagrams with a single 
insertion of this operator below. 
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limit) and it has no continuum analogue. Only when combined with ultraviolet divergences, 
(X"^CALe)ipi can acquire non-zero value. In these aspects, computation of (x^CALe)ipi 
is similar to that of quantum anomalies, although this breaking of supersymmetry is not a 
genuine anomaly in a conventional sense. ^ 

We expand F according to a number of internal loops of IPI diagrams: 

r = To + A + r2 + • • • . (3.10) 

The loop counting parameter in the present model is g*g. The tree-level effective action /q 
is nothing but the total action (3.2), 

To = -Stot. (3.11) 

In fact, it is easy to see that 5tot. satisfies the lattice Ward-Takahashi identity. In this tree 
level approximation, the breaking term vanishes in the a — > limit, 

lim x^CALe = 0, (3.12) 

a— +0 

because the Leibniz rule holds in this limit. Thus the last two lines of the identity (3.7) 
vanish in the a — ^ limit and the supersymmetry is restored in this classical continuum 
limit. 

In loop diagrams, all external lines arc tilded variables. The Ward-Takahashi identity 
for the effective action Fn (n > 1) can thus be written as 

- rn^{V2P+{Di^ + F)e + V2P_{Di^* + F*)e} 
+ E rn-tV2e^CD,P^x + E rn^^V2e^CD,P.x 



+ + L.,,e,ei;^/n + •••}= 0. (3.13) 



We also recall that our system S possesses two global U(l) "symmetries": U(l) sym- 
metry, eqs. (2.21) and (2.22), and U(l)/j symmetry, eqs. (2.24) and (2.25). In terms of the 
IPI effective action Fn {n > 1), these invariance can be expressed as 

y |r„i(^ - Fr,^6* - rn-^F + Fr,^F*\ - F^-^g + Fn^g* = 0, (3.14) 
d(l) d(f)* dF dF* J dg^ ""dg*^ ^ ' 



®It will be removed by local counter terms (supersymmetry is thought to be anomaly-free at least in 
perturbation theory) and also the structure of {x^CAL€)ipi is not universal, i.e., it will quite depend on a 
lattice formulation one adopts. 
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and 



d d ~ d ~ d ~ d 

d(P d(P* dF dF* 



-2r„^m + 2r„-^m* = 0, (3.15) 
om om* 

where we have set K = L = • • • = 0. These identities are referred to as Ward-Takahashi 
identities associated to U(l) symmetries. 

3.2 Supersymmetric Ward-Takahashi identity: Improvement and renormaliza- 
tion 

Supersymmetry is not exact in the present lattice formulation of the Wess-Zumino model. 
To achieve a sTipcrsymmetric continuum limit, wc have to apply appropriate improvement 
and renormalization to the lattice action. In particular, in the continuum limit, a renor- 
malized effective action must obey the supersymmetric Ward-Takahashi identity that is 
defined by 

- j d^x r„^{ (7^a^0 + F)e + ^^2P_ ij^d^4>* + F*)e} 

+ / d^xrn-tV2e'^CP+x+ [ d^xr„4-\/2e^CP_x 
J oq) J 0(p* 

+ J d^xrn^V2e^Cji,di,P+x + J d\rn-^V2e^Cj^d^P.x = 0, (3.16) 



and 



10 



/ 



Also the U(l) symmetries must be preserved 



I 



and 



I 



F - 4r„, 



ax 50* dF dF* 



-2r„— m + 2r„— ^m* = 0. (3.19) 
om om* 



In these expressions, all fields, external sources, mass parameters and coupling constants 
are regarded as renormalized quantities. To achieve this supersymmetric finite continuum 
theory, we take following steps. (1) We compute IPI nth loop diagrams by using the 



'^"This requirement may seem somewhat ad hoc. However, without this kind of additional condition, 

divergences arise from sub-diagrams in a IPI diagram which contains source terms K, L, M, . . . cannot be 
cancelled by an addition of local terms to Fn+i- Note that this requirement is satisfied in the tree level. 
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total action S'tot.- (2) We add appropriate local counter terms to the total action 5tot. so 
that the super symmetric Ward-Takahashi identities (3.16) and (3.17) hold in the a ^ 
limit. At this stage, all fields and parameters in eqs. (3.16) and (3.17) are understood as 
bare quantities. This improvement step removes supersymmetry breaking due to lattice 
artifacts in our formulation. (3) Then we further add appropriate local counter terms to 
the total action S'tot. so that is finite in the a — ^ limit. This step corresponds to 
the standard supersymmetric renormalization and we assume its validity, i.e., we assume 
a renormalizability of this lattice model. Explicitly, we modify the total action as 

Stot. ^ Stot. +a^Yl Zn[\x^CDx + + F*F 

X 

+FD2(t) + F*D2<j>* - -X'^CX - -(jr$ + jr*$*)|, (3.20) 

where is a common wave function renormalization factor. (4) All these steps must be 
consistent with global U(l) symmetries. We then repeat the above steps for IPI diagrams 
of one loop higher, Fn+i-, by using S'tot. so determined. 

The tree-level effective action Fq is given by the total action in eq. (3.2) and of course it 
does not need any renormalization and improvement. To see the situation in the one-loop 
level, we evaluate in the next subsection the supersymmetry breaking term in the lattice 
Ward-Takahashi identity (3.13) to this order. 

3.3 One-loop evaluation of the breaking term 

We evaluate the supersymmetry breaking term in eq. (3.13) in the one-loop order. We 
set if = L = • • • = 0. It is given by one-loop IPI diagrams with a single insertion of the 
operator (3.9). A computation of the continuum limit of these one-loop IPI diagrams is 
not difficult, if one invokes a powerful Reisz's theorem [31, 32] on lattice Feynman integrals. 
Most singular one-loop diagrams which possibly contribute to 

limV-I^Ae^ = lim|-a^y(x^CAL);^7Pe |, (3.21) 

a^O^ dK„ K=L=-=0 a^oi '^^^ K=L=-=oi ^ ' 



X 



are given by figures 1-4 and their complex conjugate. 

By using propagators in eq. (2.35), the contribution from figure 1 to eq. (3.21) is given 

by 

+2V2/5a^^a4 5;X^(x)CP+F*(y) / ^ e^''^ / ^ Xi(A;, g; m, a)6, (3.22) 

where B denotes the Brillouin zone 

B={v^R^\ Ip^I < 7r/a} , (3.23) 

and the integrand is given by 

Xi(/c, g; m, a) 

1 5i(A; + g)-5i(A;)-5i(g) 



\D2{k) + (l - l^D^ik)^ m*m f D2(A; + q) + (l - ^aD2(A: + g)) m*m 



(3.24) 
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Figure 1: A diagram which contributes to eq. (3.21). The blob indicates the supersymmetry 
breaking operator (3.9). The broken Hne is the propagator of (p and the arrow indicates a flow of 
the chirality. 




Figure 2: A diagram which possibly contributes to eq. (3.21). The bold line is the propagator 
of X- The a — > limit of this diagram turns out to be vanishing. 




I 



T 

I 

4> 

Figure 3: A diagram which possibly contributes to eq. (3.21). The a — > limit of this diagram 
turns out to be vanishing due to the U{1)r symmetry of the m = case. 

In this expression, Di and D2 denote the momentum representation of difference operators 
in eq. (2.10) 

Die'P'' = 5l(p)e*^'^ D2e'P'-' = ^2^)6*^'^, (3.25) 



-15- 



Figure 4: A diagram which possibly contributes to eq. (3.21). The a — > hmit of this diagram 
turns out to be vanishing due to the U{1)r symmetry of the m = case. 



and the explicit forms are given by 

D,{p) = l{l - (l - la^T.Pl){^ + r'i:plpl}"% (3-26) 

/i i><p 

with abbreviations ^ ^ 

Pp. = - sm{apf,), Pp, = - sin [^) ■ (3.27) 

We note that both and reduce to the momentum in the continuum theory in the 

a ^ limit, lima^oP/i = lima^o^M = Pm- 

Now, a crucial idea of the Reisz power counting theorem is to consider the A ^ oo 
limit of the integrand, after replacing the internal loop momenta by \ki and the lattice 
spacing a by a/A. From the above explicit form of the integrand, we find 

Tx{\k, q- m, a/ A) = 0(1/A^), (3.28) 

in the X ^ oo limit. This implies that the degree of divergence of the above loop integral 
(in a sense of the Reisz power counting theorem) is and the a — ^ limit of the loop 
integral may exhibit a logarithmic divergence of the form log a. 

To reduce the degree of divergence, we thus take the first term in the Taylor expansion 
of the integrand with respect to the external momentum q and consider a subtraction of 
the form 

Ii{k,q;m,a) —li{k,0;m,a). (3.29) 

However, since li{k,0;m,a) = 0, this subtraction does not improve the convergence be- 
havior at all. 

We are thus lead to consider a subtraction to the next order term in the Taylor ex- 
pansion 

Ii(k, q; m, a) — Ii(k, 0; m, a) — q^d^Ii{k, 0; m, a). (3.30) 



Then we find that this combination behaves as 0(1/A^) in the A ^ oo hmit defined above. 
The Reisz power counting theorem then states that the a of the loop integral of this 
combination is convergent and moreover the limit is given by^^ 

^^^n / 7TT4 {^i{k,q;m,a) -li{k,0;m,a) - q^dp:i{k,0;m,a)} 

= 0. (3.31) 

The last equality follows from a property of the operator AL (3.9) that it vanishes in the 
classical continuum limit. 
In this way, we obtain 

f d'^k 

= ^^^n / T^T^ {liik,0;m,a) + q^d1,Iiik,0;m,a)} 
f d^fc 

^ (2^ l™o ^^''^M^i^^/^' 0; 
= -^-^ilnQn, (3.32) 
where the coefficient ri is given by an integral 



(3-33) 

with abbreviations 

A:^ = sinfe^, fe^ = 2sin(^), 5 = l + i^gfe2_ (334) 

Prom a numerical integration, we have 

ri = +0.1518. (3.35) 
Eq. (3.32) then implies that figure 1 gives rise to 

X 

= -n9*9 J d^x [F*V2e^C-f^d^P+x + FV2e^Cji,d^P^xy (3.36) 



'^In ref. [3], this fact is referred to as an empirical "rule". 
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Next, we consider the contribution of figure 2 wliich is given by 

+4V25*5a' E E x''{x)CP^4>{y) ^ ^ e^'^^^"^) ^ ^ X2{K m, a)6, (3.37) 



where 

l2{k,q;m,a) 



Diik) Di{k + q)-Di{k)-Di{q) 

(fc) + (l - 5«52(A;)) m*m f 52(A; + g) + (^1 - ^0^2 (A; + g)) m*m 



. (3.38) 



We find that X2(AA;, g; m, a/A) = 0(1/A^) in the A — 00 hmit and the degree of divergence 
of the above integral is 1. A twice subtraction 

l2{k,q;m,a) — l2{k,0;m,a) — q^d'^l2{k,0;m,a), (3.39) 

then makes the degree of divergence —1 and the a — > hmit of the loop integral convergent. 
Again, as before, the a ^ limit of the combination (3.39) vanishes. However, this time, 
l2{k,0;m,a) = and the remaining lattice integral vanishes too: 

r (|4^ 

Jts (2^ 0; m, a) = 0. (3.40) 

Thus figure 2 has no contribution in the a ^ limit. 
The contribution of figure 3 is given by 

f d'^k 

^ / 77r~^^3{k,qi,q2;rn,a)e, (3.41) 



where 

I3ik,qi,q2]m,a) 



Di{k-q2) l-1D2{k) 



lD2{k - q2) + (1 - \aD2{k - 92)) m*m lD2{k) + (l - \aD2{k)) 



mm 



Di{k + qi)-D,{k)-Diiq,) ^^^^^^ 



lD2{k + qi) + (1 - ^aD2{k + gi)) m*m 

This behaves as 0(1/A^) in the A — > 00 limit defined above (the degree of divergence 
is —1) and the a ^ limit of the loop integral converges without any subtraction. Since 
lima^o2^3(^) Qi o-) = 0, however, the contribution of figure 3 vanishes in the a ^ limit. 
An underlying reason for this good convergence behavior of figure 3 is the chiral U(1)r 
symmetry of S with m = under eq. (2.23). Due to this symmetry, this diagram vanishes 
for m = even with a 7^ 0. In fact, if this diagram had a contribution to the breaking 
term when m = 0, through the Ward-Takahashi identity (3.13), there must be terms of the 
form (^*(/)^ or F*(p'^ in the one- loop effective action Fi: Both are however forbidden by the 
chiral U(1)r symmetry (2.24). Similarly we find that figure 4 has no contribution in the 
a — > limit. 
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As is also clear from above expressions, the ver- 
/ tex for the breaking term x^ClS.Le in cq. (3.9) in 

the momentum space (see figure 5) is proportional 
to a combination -Di(pi) + -Di(p2) + Di{p^). Prom 

^ the fact that in one-loop diagrams one of three legs 

of the vertex is always an external line and from the 
^ ^ momentum conservation at the vertex, we find that 

^ „ this vertex, when inserted in a one-loop diagram, 
gives rise to O(A^) factor in the A — oo limit de- 
Figure 5: The supersymmetry breaking fined above. This implies that this vertex effectively 
vertex (3.9) in the momentum space. behaves as a non-derivative coupling in the power 

counting argument for one-loop diagrams. Then it 
is easy to confirm that the degree of divergence of a loop integral in all possible one-loop 
diagrams other than those in figures 1-4 is negative. By repeating a similar argument as 
above, we then infer that the a — >^ limit of those contributions is zero. 

In summary, only the contribution of figure 1 survives in the a — ^ limit and we have 

lim A er> I 

a-^O^ dKa \k=L=-=0 

= -ng*g J d^x [F*V2e^Cj^d^P+x + FV2e^C^^d^P.x} ■ (3-43) 

As a general argument shows, the supersymmetry breaking in the a — ^ limit is a local 
polynomial of fields. 

3.4 One-loop level improvement and the renormalization 

Prom the above one-loop calculation (3.43), we can extract following information. Pirst, 
combining it with the a — limit of the lattice Ward-Takahashi identity (3.13), we find 
that the coefficient of the term f d^xF*F in Fi is different from the supersymmetric 
value. Namely, the coefficient of this term is not consistent with the supersymmetric 
Ward-Takahashi identity (3.16). Next, eq. (3.43) shows that Fi contains a finite term of 
the form J d'^x {F*y/2K'^C'y^dij,P-^-X + c.c.} which is not consistent with eq. (3.17). These 
two are only places in Fi in which the breaking of supersymmetry appears in the continuum 
limit (to the order 0{K, . . .)). According to the general strategy of section 3.2, we thus 
modify the total action S'tot. to 

^tot. -a^Yl n9*g{F*F + F*V2K'^CD^P^x + fV2K'^CD^P_x] 

X 

+0{K\L,...). (3.44) 

to restore the supersymmetric Ward-Takahashi identity (3.16). The added terms contribute 
to A. 

In rcf. [11], one-loop IPI two point functions are computed and it was found that 
radiative corrections to kinetic terms of x, (f) and F arc in general different by finite amount, 
although logarithmically diverging parts are of an identical magnitude. Our observation 
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above is consistent with this result and is slightly stronger: We observed that only the wave 
function renormalization of the auxiliary field F differs from other two in the continuum 
limit. The improvement above adjusts this discrepancy in wave function renormalization 
factors. 

After the above improvement, the effective action Fi becomes supersymmetric. The 

standard statements concerning the supersymmetric Wess-Zumino model are then applied 
to the one- loop effective action Fi. For example, wave function renormalization factors for 
X, <t> and F are common as indicated in eq. (3.20). Also, a local term of the form 

j d^xV{4>*A), (3.45) 

where V is an arbitrary local polynomial of and 0* without any derivatives, does not 
appear in the one- loop effective action r'l, simply because such a combination is not super- 
symmetric, i.e, it is not a solution to the supersymmetric Ward-Takahashi identity (3.16).^^ 
This conclusion is again consistent with a one- loop analysis of ref. [11] that no terms con- 
sisting only of (f) and </>* are generated by one- loop radiative corrections.^^ 

By a similar reasoning, we can also show the non-renormalization theorem [34]-[37] 
of the form in ref. [37] . The non-renormalization theorem states that the F term of the 
structure 

j d\[-x^CW"{^)P+x + FW'{4>)], (3.46) 

where W{(j)) is the superpotential, is not generated by radiative corrections, although this is 
a supersymmetric combination. We first note that, from a structure of one-loop diagrams, 
dependences of such a term on the coupling constant must be of the form 

j d^x{gFW'{g4>) + ■■■}. (3.47) 

The complex conjugate g* cannot appear here from a structure of interaction vertices. Next 
we recall that our lattice action S possesses an exact U(l) symmetry under eqs. (2.21) 
and (2.22). The above allowed structure (3.47) is, however, inconsistent with this U(l) 
symmetry. Thus we conclude that the F term (3.46) as the whole cannot be generated 
by one-loop radiative corrections. This conclusion is again consistent with the analysis of 
ref. [11]; there it was observed that terms such as x^Cx and x^C4>P+X not generated 
in the one-loop order. 

Finally, by further adding local counter terms (3.20), the one-loop effective action i^i is 
made finite, i.e., a supersymmetric renormalized theory is defined. Obviously the theory so 
defined preserves the U(l) symmetries, i.e., eqs. (3.18) and (3.19) hold, because all stages of 
the above procedure (and the lattice regular ization itself) do not affect these symmetries. 

^^An exception in this argument is V — const, that is nothing but the cosmological term. One cannot 
exclude the cosmological term from the supersymmetric Ward-Takahashi identity alone. In the present 
model, the cosmological term vanishes in the one- loop level (see appendix A), as expected from an exact 
supersymmetry of the free action So. 

^^The result of ref. [11] is somewhat stronger: Up to the quartic order in 4> or 4>* , it was observed that 
such terms are not generated even for o ^ 0. 

^^Our argument here is somewhat similar to that of ref. [38]. The conclusion here, however, is just a 
reflection of a simple fact that there is no IPI one-loop diagram made out from only F and (p external lines. 



-20- 



3.5 Higher loops 

In the one-loop order, we have observed that if we add local terms of the form (3.44) to 
the total action, especially by an adjustment of the term F*F, supersymmetry is restored 
in the continuum limit. In this subsection, we consider if this simple situation persists to 
higher orders of perturbation theory. We will find that, instead of a single combination 
of local terms (3.44), there are at most 9 combinations that we have to take into account 
for an improvement of higher loop effective action Fn- The improvement in higher loops is 
thus much more involved. We note, however, this number of combinations which require 
adjustments is much less than that in the formulation based on the Wilson fermion [3] due 
to the exact U(l)i? symmetry (2.24) and (2.25) in the present formulation. 

Suppose that the procedure in section 3.2 of the renormalization and the improvement 
(which may require an addition of 9 combinations of local terms to Stot.) work for F\, 
. .. , Fji^i. Now take a IPI n loop diagram 7„ which contains a single insertion of the 
operator (3.9). Prom the above assumption, all IPI sub-diagrams are already made finite 
by the renormalization of i^i, . . . , Fn-i- By applying a standard power counting argument 
to the present case, the superficial degree of divergence of such a diagram 7„ is given by 

<^hn) = ^-^E^-E^- 2Ef - I^F, (3.48) 

where E^, E^j,, Ep denote a number of external lines of x, ^ or (^*, F or F*, respectively 
and I(j)F denotes a number of (f)F and (i)*F* type internal lines. To derive this formula, one 
has to note that the propagator [FF*) behaves as 0(A°) (no suppression factor) in the 
Reisz power counting rule. 

Also, in deriving the above formula, we have noted a fact that an insertion of the 
vertex (3.9) in higher loop diagrams effectively behaves as 0{\) in the Reisz power counting 
rule, because in higher loop diagrams all momenta in figure 5 can simultaneously become 
large. In one-loop diagrams, on the other hand, the vertex figure 5 behaves as 0(A°) because 
one of three legs must always be an external line and momentum from the external line 
is kept fixed in the \ ^ co limit. (For one-loop diagrams, 9/2 in the formula (3.48) is 
changed to 7/2; see below.) This difference in a behavior of the supersymmetry breaking 
term in the \ ^ oo limit in one-loop and higher-loop diagrams is crucial and due to this 
difference much more combinations have to be included in local counter terms to obtain a 
supersymmetric continuum limit. 

Now, since we have assumed that all IPI sub-diagrams of 7„ are made finite by a 
renormalization, the above superficial degree of divergence will be an overall degree of 
divergence. Then if oj{'^n) < 0, the Reisz theorem states that the a ^ limit of the 
diagram 7^ is given by integrations of the a — limit of the integrand, as we have 
seen in eq. (3.31). However, due to the vertex (3.9), the a ^ Q limit of the integrand 
always vanishes. Hence diagrams which can contribute to the supersymmetry breaking in 
the continuum limit must possess <^(7n) > 0. Noting for 7„, E^ = 1, 3, . . . , we can see that 
there are seven combinations for a;(7„) > 0, i.e., {E^, E^, Ep) = (1,0,0), (1,1,0), (1,2,0), 
(1, 3, 0), (1, 0, 1), (1, 1, 1) and (3, 0, 0). The total mass dimension of 7„ is 9/2. 
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A structure of 7„ moreover must be consistent with exact global symmetries in the 
present formulation; the U(1)r symmetry, eqs. (2.24) and (2.25), and the U(l) symmetry, 
eqs. (2.21) and (2.22). The \J{1)r and the U(l) charges of are (-3, +1) or (+3, -1). 

Finally, we have to check a consistency of 7„ with the lattice Ward-Takahashi iden- 
tity (3.13). By examining various possible local terms in Fn which also must be consistent 
with the global symmetries, we finally find that a most general form of 

d 



limy^ 

a-*o ^ dKr 



K=L= 



(3.49) 



is given by 
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1 m* 

+ ^Vn— +Wn 



.m 7* 



2- 

L 5* 



• + 



m 7 



*2 



+ZnimFV2e^CP+x + m4>V2e^Cjf,df,P+x)] 



(3.50) 



plus its complex conjugate (the projection operator is replaced by P+ P-)- The real 
coefficients r„, . . . , Zn arc given by dimensionless polynomials of log{a?m*m) and a?"m*m. 

Through the a — >^ limit of the lattice Ward-Takahashi identity, the above form of the 
breaking term implies that the supersymmetry breaking terms in the effective action take 
the form: 



rr,.F*F 



+Sn4>*d J) + -i^tn4>*4) + Un' 



1*2 72 



H nVr. 



m • m 
9 9 



9 



+yn{9F4>^ + 9*F*4>*^) + Zn{mF4> + m*F* 
+{9*9r J {rnF*V2K'^Cj^d^P+x 



)} 



m 79 m 



1*2 



+ 



+ 



9 ^n*i 



A*2 ' 



1 m* 

-I o'^n— + Wn 



a? g 



m 



m 



2— 

9 9 



•*2 



m 



*2 



9 

+yn{2gF^V2K^CP+x + g4>^V2K^Cj^df,P+x) 



V2K^CP+x 



+ZnimFV2K^CP+x + m4>V2K^C-/^d^P+x) + cc.} 

+0{K\L,...). 



(3.51) 



^^To list up possible combinations, it is easier to work out the massless case m — first and then restore 
possible dependences on m by substituting 4> — > m/g and 4>* —^ m* jg* in arbitrary ways. 
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It is easy to verify that these are, modulo super symmetric combinations, most general 
local terms whose mass dimension < 4 that are consistent with the U(l)ij and U(l) global 
symmetries (the U(l) charges of the external source P±K is (^3, ±1)). This is an expected 
result from an experience in continuum theory but is not entirely trivial, because the term 
which breaks supersymmetry due to a violation of the Leibniz rule (3.9) is peculiar to 
lattice theory and cannot be treated in a framework of continuum theory. 

To remedy the above breaking of supersymmetry, we subtract eq. (3.51) from iStot. 
after transcribing it as local terms in lattice theory by substitutions 



This is the improvement step; we have to add 9 combinations of local terms to S'tot. for 
Fn to have a supersymmetric continuum limit. Then a further supersymmetric renormal- 
ization (3.20) will make Fn finite. Although this procedure may be applied in principle, 
the number of required local terms for the improvement is too many for any practical 
application of the present model. 

We can understand why the situation in the one-loop level was so simple by considering 
the case in which the first term in eq. (3.48) is 7/2 instated of 9/2. Repeating a similar 
analysis as above, as a possible form of luHa-^oY^x ai^Aea|ii"=L=-=0) we obtain 



plus the complex conjugate. This has a much simpler structure than eq. (3.50). We in fact 
found the term with r\ in the explicit one- loop calculation (3.43). We also observe that 
a general argument does not prohibit a non-zero s\; thus si = in the explicit one- loop 
calculation (3.43) is accidental. 

3.6 Exact non-linear fermionic symmetry 

Finally, we comment on an exact non-linear fermionic symmetry in this system which 
corresponds to the symmetry studied in ref. [21]. We note that the lattice Ward-Takahashi 
identity (3.6) (after setting K = L = M = -- - = Qi) may be written as 




X 





(3.52) 




(3.53) 




X 




X 



(3.54) 



where 
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However, noting the Schwinger-Dyson equation 

S-^6ip + ^ J^Jy?) = 0, (3.56) 

X X 

where ip and 6ip represent a generic field and its arbitrary variation, the identity (3.54) can 
be cast into the form 

X 

+JxS',X + J$(5^$ + J$*(5^$* + JjrS'^T + Jr*K^*}) = 0, (3.57) 

where 5'^ is a transformation modified by amount of R and TZ: 

S',x = -V2P+{Di(f) + F)e - V2P_ (£>!</.* + F*)e + Re, 

= - V2P+(Di$ + J^)e - V2P_(£»i$* + J^*)e + 7^e, 
5'^ = Se on other fields. (3.58) 

Obviously, the above form of the identity (3.57) can be regarded as a Ward-Takahashi 
identity associated to an exact symmetry 6'^ of the action S. The transformation 6'^ when 
acting on tilded variables (i.e., in the context of the Fujikawa-Ishibashi formulation) is 
nothing but the exact non-linear transformation studied in ref. [21]. However, as we have 
demonstrated, Ward-Takahashi identities in both pictures, one is based on (eq. (3.54)) 
and another is based on S'^ (eq. (3.57)), have identical contents. The presence of this exact 
symmetry itself does not imply the lattice formulation is "better" in any sense. Despite 
this exact symmetry, an adjustment of parameters is needed to obtain a supersymmetric 
continuum limit, as we have discussed in this paper. 



4. Conclusion 

In this paper, we formulated a lattice model for the N = 1 supersymmetric Wess-Zumino 
model in 4 dimensions. The U(l)ij symmetry is manifest even on a lattice with a use 
of Ginsparg- Wilson fermions. Although our formulation is perturbatively equivalent to 
the Fujikawa-Ishibashi formulation, we could avoid a singular nature of the latter formu- 
lation by introducing an auxiliary chiral supermultiplct on a lattice. Wc also analyzed 
radiative breaking of the supersymmetric Ward-Takahashi identity. The situation in the 
one-loop order is rather simple while the improvement through higher orders will be much 
more involved due to a peculiarity of the supersymmetry breaking term. In particular, in 
higher orders, we cannot avoid an adjustment of mass parameters of scalar fields which 
are quadratically diverging. In this aspect, the situation is not quite better than for the 
formulation based on the Wilson fermion. Clearly, a much clever idea is needed to achieve 
a lattice formulation of the = 1 Wess-Zumino model which avoids this too much adjust- 
ment. 

Actually, this is the starting point of the argument of ref. [7]. 
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A. Cosmological term at the one-loop level 



The cosmological term (the vacuum energy) at the one-loop level is given by a logarithm 
of the partition function with g = 0. For a correct counting of degrees of freedom, it is 
helpful to introduce real bosonic variables by 



(f>^ iA + iB)/V2, F^{F-iG)/V2, 
$ ^ (a + i(3)/V2, J^^iJ^- ig)/V2. 



Then the free part of action is represented as 



So = ^(A B F G a(5 T g)MB 



B 
F 
G 
a 

fi 
T 



X 

where matrices Mb and Mp are given by 
Mb 



fP+X\ 

P-X 
P+X 
\P-Xj 



(A.l) 



(A.2) 
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R 




R 


-/ 








--a+R 


-/ 
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(A.3) 
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where R = Re m and / = Im m and 



Mf 



2 



/D2 + m 




m 





\ 


Di 


D2 + m* 





m 




m 





-l + m 







\ 


m 





-- + m 


J 



(A.4) 



By noting the relation (2.11), we then find that gaussian integrations over bosonic variables 
give rise to the following contribution to the partition function 



^ ^^L>2 + (1- ^a-D2)m*m]}. 



-a ' 2 

On the other hand, integrations over fermionic variables give rise to 



detCdet\-^\-D2 + (1 - l:aD2)m*m\ \, 
1 4a La 2 J J 



(A.5) 



(A.6) 



where we have taken into account the fact that x and X are 4 component spinors. These 
contributions from bosonic and fermionic variables are cancelled to each other, leaving 
a constant which may be normalized to unity. Therefore, the cosmological term is not 
generated by one-loop radiative corrections even for a 7^ 0, as expected from the exact 
super symmetry of the free action Sq. 
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